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Abstract 



We obtain explicit expressions for the long range correlations in 

the ABC model and in diffusive models conditioned to produce an 

atypical current of particles. In both cases, the two-point correlation 

'^ I functions allow to detect the occurrence of a phase transition as they 

^ • become singular when the system approaches the transition. 

O 



1 Introduction 

A generic property of non-equilibrium systems, maintained in a steady state 
cn ■ by contact with several reservoirs of particles at unequal chemical potentials 

^ (or by contact with several heat baths at unequal temperatures), is the pres- 

ence of long range correlations [H El El IH El [6] . In diffusive systems, these 
00 ! long range correlations and their dependence on the system size [H [7j can be 

^ [ related to the non-local nature of the large deviation function of the density 

[8]. Another generic property of non-equilibrium systems is the possibility of 
S^ . phase transitions for one dimensional systems with short range interactions 

^: 
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[Sl[ini[IIl[I21[I31[Il[ia[ISl[ISl[ia[2D], in contrast with equilibrium sys- 
tems where such phase transitions are excluded by the well known Landau 
argument. 

The goal of the present paper is to calculate these long range correlations, 
for two examples of diffusive systems, and to show that they become singular 
at the phase transition, even when the analytic expression of quantities such 
as the average current or its large deviation function has no singularity as 
one approaches the phase transition. 

The paper is organized as follows: in section 2, we analyze the ABC model 
[T2l [T3| [m dni [18] , for which we compute first the correlation functions in the 
case of equal densities of the three species, from the known large deviation 
function of the density. We then obtain the correlation functions in the 
general case, i.e. for arbitrary densities, by two methods: a simple truncation 
procedure based on the expected scaling on the system size of the various 
correlation functions and the macroscopic fluctuation theory which was used 
recently by Bertini, De Sole, Gabrielli, Jona-Lasinio, Landim pT], [22| [23| [23] 
in the context of non-equilibrium systems. 

In section 3 we consider the weakly asymmetric exclusion process, for 
which a phase transition is expected to occur, in the large deviation function 
of the current. We calculate the pair correlation functions, conditioned on a 
current deviation, and show that they become singular at the transition. 

2 The ABC model 

The ABC model [El [131 [H [151 [H] describes a system of three species A, B, C 
of particles on a one dimensional ring of L lattice sites. Each site is occupied 
by one and only one of these three types of particles. The dynamics are fully 
specified by the following exchange rates between neighboring sites: 



(1) 



where q < 1 (the range g > 1 is related to the case g < 1 by a left-right 
symmetry). Since the number of particles of each species is conserved by 
the dynamics, the properties of the system depend on the global densities 
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fA^fB^ fc of the three species and on the parameter q. As all sites are occu- 
pied by one of the three species, one has of course 

TA + TB + rc = 1 ■ (2) 

When g = 1 all allowed configurations are equally likely (and the steady state 
density profiles of the three species are fiat) but as q decreases, the different 
species tend to segregate and in the limit g ^ the steady state consists of 
three clusters: a cluster of all A's followed by a cluster of all i?'s which is 
itself followed by a cluster of all C's. 

Let us first briefly recall some known properties of the ABC model [H]. 
As q varies, with the following scaling depending on the size L of the ring. 



q = exp 



L 



(3) 



one observes (for large L) phase transitions, in the steady state, at some 
critical value Pc- 

For equal densities (r^ = r^ = re = |), it is known that 

Pc = 27rv^ (4) 

and that for j3 < j3c the density profiles of the three species are fiat, while 
they become space-dependent for /3 > Pc- 

A stability analysis of the fiat phase leads to the following prediction of 
Pc in the general case 

Pc = T77 (5) 

[l-2{rl + rl + rl)f' 

which should be the exact expression as long as the transition remains a 
second order phase transition. One expects, however, a first order phase 
transition to occur for some range of densities at least when r\ + r\ + r"^ < 
2(r^ + r'% + r^), and the precise location of this first order transition is not 
known. In the fiat phase, the steady state currents have expressions 

Ja = P ^ 

(and similar expressions for Jb and Jc) which do not show any singularity 

as P ^ P~. 



In the steady state, the probabihty of observing macroscopic profiles 
Pa{x), Pb{x), pc{x) (where < x < 1 is a macroscopic coordinate along 
the ring) has the following large L dependence 

Fto{Pa{x), Pb{x), pc{x)) ^ exp (^- LJ='[pa{x), pb{x), pc{x)]j . (6) 

The large deviation function JF of the density profiles is known for the ABC 
model when r^ = tb = re = ■^: 

T[pAix),pB{x),pc{x)]= (7) 

K+ dx [pa (x) In Pa (x) + pb (x) In pb (x) + pc (x) In pc (x)] 
Jo 

+P dx dz z[pb{x)pc{x + z) + pc{x)pa{x + z) + pa{x)pb{x + z)] . 
Jo Jo 

where k is a normalisation constant but it is not known for arbitrary rAjTBjrc- 
So far it has only been computed in [T3| to order /?^ 

^[pAix),pBix),pcix)] = 

K+ dx[pA{x)\npA{x) + pB{x)\npB{x) + pc{x)\npc{x)] 
Jo 

+P dx dz z [pb{x)pc{x + z) + pc(x)pA(x + z) + pa{x)pb{x + z)] 
Jo Jo 

--/3M dx j dz z{l- z)[rA{l-^rA)pB{x)pc{x + z) 
4 JO JO 

+rij(l - 'irB)pc{x)pA{x + z)+ rc(l - ^rc)pA{x)pB{x + z)] 

+0(/3^) (8) 

2.1 The correlation function in the equal density case 

When the large deviation function is known, as in (jTj) for the equal density 
case, the calculation of the correlation functions in the flat phase can be 
easily done from the contribution of the fluctuations of density around the 
flat profiles. For a small perturbation of the form 

Pa{x) = — h Oi cos(27™x) + a2 sin(27rna;) 

Pb{x) = - + bi cos{2'7inx) + 62 sin(27™a;) 
o 



with n > 1 and pc{x) = 1 — Pa{x) — Pb{x), one gets from ([7]) at quadratic 
order in Oi, 02, 6i, 62 



Pro (pa (a;), Pb(x),pc (a;)) 
exp 



^-^^ 2 , ^ I, , ^2 , ^2 , ^ t I ^2n , 3I//3, 



— -(a^ + ai^i + 6^ + 02 + 02^2 + ^2) + 1 — (0'1^2 - 02^1 
/ 4n7r 

Therefore, the covariances in the steady state of the n-th Fourier coefficients 
of the densities Pa{x) and Pb{x) are 

(a?) = (al) = {h\) = {hi) ^^"'"' 



3L(12n27r2 _ ^2) 
(0102) = (6162) = 
(ai6i) = (0262) = 



^n^Tx'^ 



(ai&2) = -(02^1) 



3L(12n27r2 _ ;52) 
47™/3 



3L(12n27r2 - /?2) 

Then if one sums over all the fluctuations, i.e. over all the wave numbers n 
on the ring, one obtains 

and similar expressions for the other correlation functions (we will see be- 
low in fl36|37l) a different way of writing these expressions). Clearly these 
expressions become singular as /? approaches the critical value Pc = 27rv^. 

2.2 The correlation function for arbitrary densities by 
a simple truncation procedure 

Let us introduce the notation Ai = 1 when site i is occupied by a A particle, 
Ai = otherwise, and similarly Bi or Cj for a i? or a C particle, so that on 
each site i one has 

A + B, + Q = l. (11) 



One can derive from the dynamical rule ([T]) the following exact evolution 
equations: 

at 

+ (A_ia) + g(aA,+i) - (A,a+i) - g(Q_iA,) , (12) 

where ( ) stands for the expectation with respect to the dynamics. For 

|J-^I>2 



djAjAj) 
dt 



q{Ai_iBiAj) + {B,A+iAj) - q{ABi+iAj) - (B^^iAA^) 

+q{AiAj^,Bj) + (AiBjAj+i) - q{AA,Bj+i) - {AB^^iA,) 
+ (A,_iC,A,) + q{CiAi+iAj) - {Aid+iAj) - q{a^iAAj) 
+ {AA,_,C,) + q{AC,A^,) - (A,A,C,+i) - q{A,C,.,A^) (13) 



djAB,) 
dt 



= q{A^iB,Bj) + {B,A+iB,) - q{AB,+^Bj) - {B^^iAB,) 

+q{AB,.,C,) + {AC,B,^i) - q{AB,C,+i) - {AC^^^B,) 
+ {A^,QB^) + q{QA^,B^) - {AQ+,B,) - q{Q^,AB^) 
+ {ABj^iAj) + q{AAjBj+i) - {AB^Aj+i) - q{AAj^iBj) (14) 

and for j = i + 1: 

'Ml^^ = q{A^,B,B.,^,) + (5,A+i) - q{A,B,+i) - {B,^,A,B,+,) 

+ {A^Ci+iB.i+2) - q{AB.i+iC.i+2) 
+ {A^iCiBi+i) - q{a^iABi+i) + q{AA+iB,+2) - {ABi+iA+2) 

(15) 

and similar expressions for the evolution of the other one or two point func- 
tions. In the steady state, the left hand sides of the equations (I12m5p vanish. 
For diffusive systems, one expects the long range part of the p-point 
connected part of the correlation functions to scale as L^~p in the steady 
state [H El E]. This means that, for large L, in the flat phase of the ABC 
model, 

{AB,) = {A,){Bj) + {ABj), = TATB + -^Fab (^) (16) 
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where Fab stands for the normahzed connected two-point function and 



{A.Bfik) = TATBrc + 



1 

L 



^ I i — i\ ^ ( k — i\ „ ( k — j 



L J ^ ^^\ L ^ "^\ L 



where Gabc stands for the normahzed connected three-point function. Using 
similar notations for ah the other 2-point and 3-point functions and the 
relation fITT]) . and replacing the various correlation functions in fll3|14lll5p by 
their scaling forms fll6lll7p . one gets to leading orders in L~^ (order L~^ for 
fITMD . and order L~^ for (|T^) 

Faa{xy = f3rA[Fba{xy - Fab{x)'] 

Fbb{x)" = l3rB[Fba{x)' - Fab{x)'] (18) 

Fabix)" = (3[rBFaa{x)' + r^F,,(x)' - (1 - 2rA - 2rB)Fab{x)'] 
Fba{x)" = l3[-rBFaa{x)' - rAFbb{x)' + (1 - 2rA - 2rB)Fba{xy] 

as well as 

Fab{Q)-Fba{Q) = ?>PrArBrc. (19) 

Note that in flTSl) at leading order in 1/L, the contribution of the functions 
G"s which appear in flTTl) drops out. Note also that 

Fab{x) = Fbail - X) . 

The general solution of (fT8l) is 



Faaix) = Kaa ' 2i^/3r^ COS (^^(Sx 



Fbbix) = Kbb - 2K(3rB cos ( ^^(2x - 1) j (20) 

Fab[x) = Fba{l -x) = Kab + K{1 - 2rc)/3cos f^(2x - 1)) 



-K(3^sm ( ^^(2x - 1) J 

where K, Kaa, Kab, Kbb are integration constants and 

A = l-2{r\ + rl + rl). (21) 



Relation fll9p leads to 

K = -^l^L^l^ . (22) 

The remaining constants can be determined by using the fact that 

k=l 

and similar sum rules for the other correlation functions which imply that 

/ Faa{x)dx = -rA{l-rA), / Fbb{x)dx = -rB(l-rB), / Fab{x)dx = TaTb 
Jo Jo Jo 

and the final result is for < x < 1 

Faa{x) = -rA{l-rA) — ^-7^ I /?cos ( -^^{2x 

V A sin ^=^ 




Fatix) = TATB + ^^ 



2VAsin 2 



2rc) /?cosp^(2x-l) -2- ^ ' ' 



-/?VAsin f ^ 



2 (2x-l) (24) 



As /5 ^ ^ the correlation functions diverge and one recovers the ex- 
pression ([5]) of the second order transition. 

As starting from ( TTGfTTj) we only considered the correlations in the flat 
phase, our truncation procedure is unable to distinguish between a stable and 
a metastable phase and to predict the occurrence of first order transitions. 
Thus we cannot exclude that the fiat phase becomes metastable for some 
values of [3 smaller than (3c given by ^ . 

2.3 The macroscopic fluctuation theory 

In the ABC model there is always one and only one particle per site. One 
can then describe the system on a macroscopic scale and on a diffusive time, 
i.e. on times which scale as L^, by only two density profiles pyi(x, t), pb{x, t). 

8 



The key quantities to the macroscopic description of the system are the 
typical currents for given density profiles and their variance in the steady 
state. Let Qf{i),Qf{i) be the fluxes of A, B particles between sites i and 
i + 1 during a long microscopic time t. Starting at time from smooth 
macroscopic profiles pa{x), pb{x), the densities evolve according to [1] 



drPA{x, t) = dlpA{x, t) + f3d^{pA{x, t){pb{x, t) - pc{x, r))) (25) 

drPB^x, r) = d^psix, t) + pd^ipsix, t){pc{x, r) - pa{x, r))) 

where the macroscopic time r scales like t/L"^. One can understand (125!) 
from (fT2l) by saying that for large L, there is a local equilibrium at position x, 
characterized by the densities Pa{x, t) and Pb{x, t). In this local equilibrium, 
the hydrodynamic currents are given by 

U = —^ - /3pa{pb - Pc) (26) 

(IB = —^ - Ppb{pc - Pa) 

To compute the covariance matrix of the currents, it is enough to consider the 
system in equilibrium with /3 = and constant densities pA-, Pb- In this case, 
the A-particles evolve as a Symmetric Simple Exclusion Process (SSEP) (the 
B, C particles play the role of holes) and the variance of the total current 
Qt^ = EiiQf(Oisgivenby[n] 

aAA = }i^^{{Qf?) = 2pA{l-pA). 

t^oo tL 

In the same way, the variance of Qf = Y.i=i Qf (0 is given by asB = 2pb(1 — 
Pb)- Finally, if we view the C-particles as holes and the species A,B as a. 
single type of particle, then the evolution is the same as the SSEP at density 
Pa + Pb, so that 

lim —{iQt + Qff) = 2{pA + Pb)(1 -PA~ Pb) = (Taa + (Tba + 2a ab ■ 

where the covariance between the currents of A and B is then aAB = —2paPb- 

One can adapt the macroscopic fluctuation theory [31], [32] to predict that 
the probability of observing density profiles pa{x, t), pb{x, t) and (rescaled) 



currents Ja{x,t),Jb{x,t) over a (rescaled) time interval < r < T is given 
(this rescaling is discussed for example in section 5 of 



where 



with 



Fto{pa,Pb,Ja,Jb) ~ exp i^-L I[o^t]{pa, PbJaJb)) , (27) 

^Io,t]{pa,Pb,Ja,Jb) = J dx j dTn{pA,pB,JA,JB), (28) 



2 



_ (^BBJJA - QaY - '^(^AbHa - qA)iJB - Qb) + (^AaJJB - Qb] 

'2'{(Taa(^bb - (^ab) 
Note that instead of fl29|) . one could write Langevin equations for the currents 



dpA ^ . ^ , Va{x,t) 

U = -^-Ppa{pb-Pc) + ^^ 

dpB n I N , r]B{x,T) 

JB = — 1 PPb[Pc - Pa) + 



dx a/L 

with the following correlations of the white noises riA{x, r) and rjBix, r) 

{r]A{x, t)7]b{x, t')) = (Jab{pa{x, t),pb{x, r)) 6{r - r') 

and similar expressions for (?7yi(x, r)?7^(x, r')) and {r]B{x,T)riB{x,T')). 

As always, the currents and the densities are related by the conservation 

laws 

dpA ^ _dJA . dpB^ ^ dJB 

dr dx ' dr dx 

For a perturbation of the form (which satisfies these conservation laws) 

Pa{x, t) = rA + k{ai cos[A;x + lot] + a2sin[A;x + cur]) 
Pb{x, r) = tb + k{hi cos [/ex + ujt] + 62 sin[fc2; + ut\) 
j^(x, r) = —[3rA{rB — re) — uj{ai cos[kx + ujt] + 02 sm[kx + ut\) 
Jb{x, r) = —PrB^rc — r^) — uj{bi cos[kx + ujt] + 62 sin[A;x + ujt]) (30) 

where 

2rmr 
k = 2mi uj = (31) 

10 



with n & N and m & Z one expands fl28H29l) to second order in ai, a2, &i, 62 
and the probabihty (1271) becomes a Gaussian which leads to the following 
expectations 

(0102) = (fel&2) = 



(«i) = («2) = -^^(/5'/cl4(l - 4r^)rBrc; - (1 - 2r^)^(l - r^)] 

+2/3fc(l - 2r^)(rc - rB)a; - (1 - r^)(A:^ + u^)) 
{a,b,) = (0262) = -^^^( - /5'A;2[4r^r5 - 1 - 2rc + Sr^] (32) 

+2/5A;(rij - taP + (A;^ + u^)) 

, , , , , , 1 24:pk^rArBrc 
(0162) = -(a2&i) = -^ ^ 

where A is defined in (12T|) and 

T, = p^k^A^ - 2p^k''A{k^ - u^) + {k^ + 0^2)2 . 

In one adds up the contributions of all perturbations of the form fl30|l . 
one gets for the correlation functions 

{pAix)pAiy))c = Y. ^^('^1) cos[k{y - x)] 

k,u) 

{pA{x)pB{y))c = J2 ^^(^1) (^os[k{y - x)] + k'^{aib2) sm[k{y - x)] 

k,uj 

In the long time limit the sums over the discrete frequencies fl3Tl) can be 
replaced by an integral and one gets 

{pAix)pA{y))c = (33) 

1 ^ rA[(3\3rAil - 2r^)2 - A(2 - 5r^)) + STr-'n^l - r^)] 

L 5, An^n^ - (3^ A ^^^^^^'^^ " "^^ 

n>l ' 

{pA{x)pB{y))c = 2 2^ 47r%2 _ ^2 A cos(27™(y - x)) 

(34) 
11 



In the particular case of equal densities, one recovers 



Using the identities 

2 ^ cos[27m(x — y)] = — 1 + 5{x — y) 



n>l 



cos[a(i — x)] cos{ax) + cos(a(l — x)) 2 -^ 



asm ■2 



asm a 



sin[a(| — x)] cos(ax) — cos(a(l — x)) y-^ 



sm 



a 



1 — COS a 



a^ ^ Air^n^ — a 
8mi 

n>l 



(35) 



cos 2mix 



Att'^u'^ — a^ 



sin 2mTX 



[1 - 2rAf + A = Arerc 



one can rewrite (13311341) as 

{pA{x)pA{y))c = 



1 

L 



rA{l - rA){S{y - x) - 1) - SPr^rsrc 



/cos(^(l-2|/ + 2x) 



V 



A sin 



PVa 



2 



(36) 



{pA{x)pBiy))c 



1 

L 



3 / COS 

-rArsiKy -x)-l)-V -PrATBrci'^ - '^rc) 



3 sin(/3^(l-2y + 2 



X 



sm ' 



^il-2y + 2x) 



PVa 



A sin 



which are totally equivalent to the expressions (I23|24l) . 

As mentioned before, the large deviation function ([6]) of the density is 
only known (jHj) at order /3^ for the ABC model. The knowledge (I36|37l) 
imposes constraints on the large deviation function. Trying to generalize ([8]), 
we found the following expression for the large deviation functional ([6D 

J='[pAix),pB{x),pc{x)] = 



2 



(37) 



(3J 
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K,+ dx[pAix) \npAix) + pb{x) hipB^x) + pc{x) In pc(a;)] 
Jo 

ri pi sin(c( — z]] 
-(3 j dx dz — .^ — [pa{x)pb{x + z) + pb{x)pc{x + z) + pc{x)pa{x + z)] 

'•1 /-i cos(c(|-^)) 



-3/3^ dx dz — - — ^— — — [rc(l - 3rc)pA{x)pB{x + z) 
Jo Jo 4csm(|) 

+rA(l - 3rA)pB{x)pcix + z) + rB{l - 3rB)pcix)pAix + z)] 
where 



c = /3y A - 9rArBrc 

was compatible both with (I8j) at order /J^ for arbitrary density fluctuations 
and with the exact expressions of the correlation functions fl36|37p . i.e. for 
arbitrary /3 but small deviations of the density from the flat profiles. What 
this expression becomes for general f3 and arbitrary deviations of the density 
is an open question. 



The macroscopic fluctuation theory developed by Bertini et al. 
(see also [25]) relates the steady state large deviation functional (^ to the dy- 
namical large deviation functional fl27H28l) . One can try to apply this general 
procedure to the ABC model. We fix < /? < /3c such that the steady state 
profiles are flat (equal to tajTb) and such that the hydrodynamic equations 
f l2S]l relax to r^, r^ for any initial density profiles (this property excludes situ- 
ations with metastable profiles). Then the density large deviation functional 
([6]) is given by 

J='[Pa{x),pb{x),pc{x)] = inf |j[„oo,o](p,j)| , (39) 

where the infimum of the functional X[_oo,o] (EH]) is taken over the trajectories 
(p(x,t),j(x,t)) = {pA{x,t),pB{x,t),]A{x,t),]B{x,t)} starting at time -oo 
from the steady state configurations {ta, tb) and ending at time t = at the 
configuration (p^(x),pb(x)) {with. pc{x) = 1—pa{x)—pb{x)). To recover the 
large deviation functional in the steady state, there is no other constraint on 
the currents j besides the conservation law dtp = —dx]- Finding the optimal 
trajectory is in general an open problem which has been solved only in rare 
instances [22| [27] . For the sake of completeness, the equations satisfied by 
the optimal trajectories are presented in Appendix II. 
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A perturbative analysis of the variational problem fl39p leads to an ex- 
pansion of T in the vicinity of the steady state. For example fl38p can be 
obtained as the second order expansion of T and the two-point correlations 
can then be recovered by inverting the quadratic form in fl38l) . Expanding 
T at the second order is essentially equivalent to the mode decomposition 
implemented in this section fl5U]) . To see this, we first note that fl5I?]) can be 
rewritten as a variational principle on the time interval [—00, 00] with the 
constraint that the trajectory p is equal to (p^(x),pB(a;)) at time 

•^[Pa (a;), Pij(x),pc' (2;)] =inf{2^[-oo,oo](p,j)} • (40) 

This identity comes from the fact that the contribution for positive times 
is null as one can choose p which relaxes to equilibrium according to fl25p . 
For small deviations, X[_oo,oo](p,j) can then be expanded with respect to the 
Fourier modes fl30l) . The different modes decouple and we recover a quadratic 
form with coefficients related to the correlations fl32l) . The constraint on the 
density at time is then achieved by optimizing this quadratic form over the 
temporal modes lo. This leads to a second order expansion of T (in terms of 
the spatial Fourier modes) from which the steady state two-point correlation 
functions (l33ll34D can be obtained. 



3 Diffusive systems on a ring 

The second example we analyze is the phase transition [291 [HI [2H] which 
is expected to occur in the large deviation function of the current in some 
diffusive systems such as the weakly asymmetric exclusion process (WASEP). 
The WASEP on a ring is defined as follows: one considers a fixed number 
N = Lr oi particles on a ring of L sites, with at most one particle per site. 
Each particle hops to its neighboring site on its right with rate q~^ and to 
its neighboring site on its left with rate q, provided that the target site is 
empty. As for the ABC model, we consider here the diffusive regime where 
q scales with the system size L as in ([3]) 



q = exp 



(1 
L 



For general diffusive systems with one species of particles on a ring, one can 
write an expression very similar to (12711291) for the probability of observing a 
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density profile p{x,t) and a (rescaled) current j{x,t) over a (rescaled) time 
interval T (the microscopic time t = TL'^) 



where 



Pro(p, j) ~ exp {-L J[o,t](p, j)) (41) 



with 



X[o,T](p,j)= / rfa;/ drHipJ), (42) 

JO JO 

^. . [j(a:,r) + D(p(a:,r))%l)-Mp(^,r))P 

and the current j and the density p are related, as usual, by the conservation 
law 

dr dx 



We refer to [28l[29] for further details on the derivation of (H2|) . The functions 
D{p) and a{p) in (H3|) are two functions characteristic of the diffusive system. 
For the WASEP, as defined above, these two functions are D{p) = 1 and 
(7(p) = 2p(l — p). In general their ratio is related to the compressibility at 
equilibrium by the Einstein relation (see [30] equation (4.6.1) or [31]). 

If one considers the flux Qt{i) of particles between sites i and z + 1 during 
a long microscopic time t, the large deviation function G{Jq) of the current 
describes the distribution of Qt{i), in the long time limit 



Pro^ = -U.exp 



-{g{jo) 



(45) 



(the L dependence which appears in ( l45l) can be understood because the 
system is diffusive [28| [29]). For systems with conservative dynamics such 
as the WASEP, which are irreducible Markov processes with a finite number 
of states, the large deviation function G{jo) depends neither on the initial 
condition at time nor on the final configuration at time t nor on the section 
i chosen. 

In what follows we consider (instead of the flux Qt{i) through a section) 
the total flux Qt 
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which has the same large deviation function 

(47) 



Pro ( Y = io ) ~ exp 



-jGUo] 



According to the macroscopic fluctuation theory, one can calculate this 
large deviation function G(jo) by looking at the time evolution of the density 
and of the current which maximizes (HTl) with the constraint that 



T rl 



SO that 



dr dx j{x,t) = jo (4J 

1 JO Jo 



GUo) = lim 7^ mill J[o,t](p, j) (49) 

J — »oo J \ p{x,t),j(x,t) I 



When the optimal density and current profiles in (l49l) are constant in 
time and in space [29], the expression of G{jo) follows immediately 

where r = N/L is the density of particles along the ring. 

By analyzing the neighborhood of the fiat profile solution, it was shown 
in [16] that this profile is locally unstable when 

87r^D\r)a{r) + [(3^a\r) - j^]a"{r) < (51) 

and the optimal profiles becomes space or space and time dependent. (Note 
that, as for the ABC model, (15T]) is obtained by a local stability analysis. 
Thus the optimal profiles might already be space or space and time dependent 
even when (1511) is not satisfied with the occurrence of first order transitions). 

We are now going to see that, although there is no trace of the second 
order phase transition flSTl) in the expression (l50l) . the correlation functions 
become singular along the stability line (|5T!) . There are several ways of 
defining the correlation functions, conditioned on the current: if the flux is 
conditioned to be Qt = jot over a time interval —t/2 < t' < t/2, one expects 
the correlation {p{x,t')p{y,t'))c to depend on whether t' is close or far from 
the boundaries of the time interval (—t/2, t/2). We are now going to calculate 
two of these correlation functions 

(p(x)p(y))intermediate = hm (p(x, 0)p{y, 0)\Qt = jot)c (52) 
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(p(x)p(y))finai = lim (p(x, t/2)p{y, t/2)\Qt = jot)c (53) 

t — ^oo 

During the macroscopic time interval [— T/2,T/2], one considers a small 
perturbation of the flat profile r of the form 

p(x, t) = r + k[fi{T) cos{kx) + /2(t) sm{kx)] (54) 

j{x, r) = /o(r) + /2(r) cos(A;x) - /((r) sin(A;x) 

where /i(r) and /2(t) are a priori two arbitrary functions and k = 2Txn with 
77, > 1, integer). The parameter Iq{t) takes into account the constraint on 
the current in the time interval [— T/2, T/2] 

r /o(r)=jo, if -T/2<r<r/2, 
1 /o(r) = /5(T(r), otherwise. 

One can rewrite (143|) to quadratic order 

7^(p, j) = A(r)(A^+/2^)+5(r)(/f +/^^)+C(r)(/i/^/2/0+i?(r)(/i/{+/2/^) 

(55) 
where 

^^ = —A + O + 7^^ O 9 

4(7 8 40-3 80-2 

5(r) = -, C(r) = - '"^y , i^(r) = ^ (56) 

^ ^ 4(t' ^ ' 2^2 ' ^ ^ 2a ^ ^ 

with the functions D, a, a', a" evaluated at p = r = N/L. 

When T goes to infinity, the expectations (/f(0)), (/i(0)/i(0)), (/|(0)) 
can then be computed in the following two situations (see Appendix I): 

• in the case fl52l) . one has /o(t) = jo for all times r. One gets from flMI) 
(/i(0)^) = (/2(0)^) = , X, „ ,,^,,.. (57) 



2 2o- 

and by summing the contributions of all the modes, this gives 

27r .^-^ n cos(27rn(y — x)) a 

{P{x)p{y)) intermediate = -T 1^ , „, „ ■ (58) 

V z zo" 
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in the case (13^ . one has Io{t) = jo for r < and Io{r) = [3a{r) for 
r > 



(/i(0)^) = (/2(0)^) 





2a 




L 


Dk^ + ^iD^k^ + ^'''Y' - 


2a 



(59) 



and the sum over all the modes leads to 

47r ,-^ ncos{27m{y — x)) a 



{p{x)p{y)) 



final 



^E 



i nil 2nnD + Jin^n^D^ + E^ - ^ 

V z za 



(60) 



We see in expressions f l58|60l) that the correlation functions become singular 
along the transition line (ISTjl . where the n = 1 mode becomes unstable. Note 
that only the intermediate correlation function flSSl) diverges at the transition. 



Let jo = Pclr) be the mean current, then one can rewrite the correlation 
function (158|) as 



{p[x)p[y)) intermediate 

a 



(61) 



2LD 



^2cos(27™(y-x)) 1 



n>l 



8cr7r2n2Z)2 



(Jo^ - (Jo)^) 



When there is no constraint on the current then jo = Jo and (!6T]) reduces 
(see (BSD) to 



(p(a;)p(y)) 



cr 



intermediate 



2LD 



y^ 2 cos{2'Kn{y — x)) 



a 



n>l 



2LD 



{5{x-y)-l), 



which are the correlations for the (micro-canonical) invariant measure of the 
diffusive part of the dynamics. This was already emphasized in [6l [25] . On 
the other hand, if the system is conditioned to an atypical current deviation, 
the term in the bracket in (!6T|) does not vanish. It is interesting to note that 
all the modes have the same sign which depends on Jq — Jq and a" . One can 
check that depending on the sign of Jq—Jq the correlations for x — y small are 
either positive or negative: the particles tend to cluster or to spread according 
to the constraint on the current. We interpret this clustering phenomenon as 
a precursor of the macroscopic clustering which occurs after the transition. 
Instead of the perturbation (l54l) . one could also have used a space/time 
decomposition over the Fourier modes as in ( l30l) . However this would have 



led to extra difficulties to treat the fluctuations at the final time. In the latter 
case, the functional depends on the current constraint which is not uniform 
in time so that the contributions of the modes at different frequencies would 
be coupled (unlike the case of intermediate fluctuations). In the particular 
case of the WASEP, a microscopic approach based on a truncation procedure 
as in section [22] would also have led to the correlations (ISSjl . 

In [in], it was shown recently that a very similar phase transition occurs 
in the SSEP (the symmetric exclusion process) when one considers the large 
deviation function of the activity Kt (which is the total number of changes 
of conflgurations during time t). A calculation almost identical to the one 
presented in this section can be done to calculate the correlation functions, 
conditioned on the value of Kt- The nature of the singularities of the corre- 
lation functions are then very similar to those found in fl58| [60l) . 

4 Conclusion 

In the present work, we have calculated, using two different approaches, the 
long range correlation functions in the ABC model and shown that they 
become singular at the second order phase transition fl9|10ll33ll34|36|37l) . We 
have also calculated, for general diffusive systems, the two point correlation 
functions, conditioned on the current flSSj 1501) and seen that they become 
singular at a second order phase transition. Similar calculations can be done, 
when conditioned on other quantities such as the activity [19]. 

It would be interesting to try to extend our results to other phases than 
the flat phase or to other situations than the ring geometry, like open systems. 
The large scale Gaussian fluctuations of the density, which here are at the 
origin of the long range correlations, allow also to calculate all the cumulants 
of the current [19] for diffusive systems. It would be interesting to see whether 
one could establish more direct or more general relations between the long 
range correlations and the distribution of the current, and to know how these 
relations are modifled in the case of open systems. 

Lastly we have obtained an improved expression fl5S]) of the large devi- 
ation functional of the density for the ABC model. This expression is still 
an approximation. It would be interesting to go further, for example by 
calculating higher correlations using the truncation procedure ( TT61IT7I) . 
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Appendix I 

In this appendix, we establish the expressions fl57|59p . Let us consider two 
Gaussian functions /ii(T), /i2(t) of the time r distributed according to the 
following distribution 



Pro({/ii(r), h2{r)}) ~ exp 



/oo 
dT[A{T){hl + hi) + B{T){h[^ + h'^) 
-oo 



+C{T){hh'^ - h'.h) + E{r){hh[ + h2h'^)] (62) 

Here we are interested in situations where the functions A{t), B{r), C{t), E{t) 
take some constant values A., i?_, C_, E^ for r < and other constant val- 
ues (possibly the same) y4+, 5+, C+, i?+ for r > 0. 

If we fix /ii(0) and /i2(0), and one integrates (162!) over /ii(t) and h2{T) for 
all times r 7^ one gets 



Pro({/ii(0),/i2(0)})~exp 
where 



771 771 

-L ( a+5+ + a_5_ + — ) (/ii(0)2 + /i2(0)^^ 



AA,B,-Cl \AA^B_-C_ 



+ 2B+ ' 2S_ 

In fact it can be shown that the functions which maximize ( l62l) for r > 
are, at fixed /ii(0), /i2(0) 

/ii(r) = e-"+"[/ii(0) cos(/5+r) + /i2(0) sin(/?+r)] (63) 

/i2(r) = e^"+^[/i2(0) cos(/3+r) - /ii(0) sin(/5+r)] 
where 



25+ 

Similar expressions with a+ and /?+ replaced by — a_ and (3- give these 
optimal functions for r < 0. This leads to 

'"■'"^ - <"^<°« - |2(a.,fl,+a_i;'_) + g_-i;,|L (''^' 

and (/ii(0)/i2(0)) = 0. 
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Appendix II 

In this Appendix, we study the variational problem (l39l) . We set 

^(f,\ _ ( 2Pa(1 - Pa) -'^PaPb \ j^^.. _ ( Pa{pa + 2pB - 1) \ 

^^^^ 1^ -2pAPB 2pb{1 -Pb) ) ^^' \ Pb{1 - 2pA - Pb) ) 

The trajectories are denoted by p = {pa, Pb) and the currents j = {Ja,Jb)- 
It will be convenient to rewrite the currents in terms of the new variables 
H = {Ha, Hb) 

j(x, t) = -dj{x, t) - PF{p{x, t)) + S(/3(x, t))d,H{x, t) , (65) 

with H{0,t) = H{l,t) = 0. In particular 

dtp{x,t)=dlp{x,t)+(3d,[F{p{x,t))-d^{^{p{x,t))d^H{x,t))). 

Thus the functional (l28l) reads 

2:[_oo,o](p,j) = \f^Jtf^d^ (H'a.K) S(p(x,t)) ( fM ' (66) 

where H'^, H'^ stand for the spatial derivatives. Optimizing the functional 
implies that the optimal trajectories satisfy 

dtPA = ^PA + (3dx(^pA{pA + 2pB - 1)) - 29x.(pa(1 - Pa)H'a - paPbH'b) 

dtPB = ApB + (3d^{pBil - 2pA - Pb)) - 29,.(pb(1 - Pb)H's - PaPbH'a) (67) 

dtHA = -AHA + P{2pA + 2pB-l)H'^-2f]pBH's-{il-2pA){K)'-2pBH'^H's 

dtHB = -AHB + 2f]pAK-(3{2pB + 2pA-l)H's-({l-2pB)iH's)'-2pAH'^H's 



with the constraint on the densities at time 0, {pa{x, 0), pb{x, 0)) = {pa{x), pb{x)) 
and {va, tb) at time — cxd. 

Solving fl671) amounts to knowing the functional T (see [22] equation 

(2.16)) 

H -^^ H -^^ 

nA — 7^ — , ns — 7^ — . 

OpA OpB 
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In the case r^ = r^ = 1/3, the optimal drifts can be guessed from the explicit 
expression (jTj) of the functional T 



(0 A {x\ \ f t^ 
TT +/^ / <iupB{u)- pa{u)^ 
I - Pa{x) - pb{x) J \Jo 

du (1 - Pb{u) - Pa{u)) + / du{u- x){3pb{u) - 1) j 



and 



TT, I N {Pb{x) - 1)p'a{x) -pa{x)p'b{x) f 

nAx) = ——. — — ^ JiSpBix) - 1 



Pa{x){pa{x) + Pb{x) - 1 
A similar expression holds for Hb leading to 

dtpA = -ApA{x,t) - (3d^(^pA{pA + 2pB - 1^ 
dtPB = -/^PB{x,t)- (3dJpB{l-2pA- Pb 



Thus the optimal trajectories are the time reversed of the hydrodynamic 
equations (l25ll . This could have been guessed from the microscopic reversibil- 
ity derived in [13]. Under the assumption that the hydrodynamic equations 
fl25|l relax to the fiat densities r-^ = tb = 1/3, the expression ([7]) of JF can be 
recovered from the optimal trajectories by integrating ( l66l) over the time. 
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